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A few recent works suggest the possibility of controlling light propagation at the interface of periodic multilayers 
supporting Bloch surface waves (BSWs), but optical resonators based on BSWs are yet to demonstrate. Here we 
discuss the feasibility of exploiting guided BSWs in a ring resonator configuration. In particular, we investigate 
the main issues related to the design of these structures, and we discuss about their limitations in terms of 
quality factors and dimensions. We believe these results might be useful for the development of a complete 
BSW-based platform for application ranging from optical sensing to the study of the light-matter interaction 
in micro and nano structures. 
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1. Introduction 

Light confinement near the surface of a photonic struc¬ 
ture is crucial in a number of applications involving the 
light-matter interaction. For instance, it is appealing 
for the realization of optical sensors, in which one re¬ 
lies on the interaction between the electromagnetic field 
and the analytes under investigation. In this respect, 
surface plasmon polaritons are probably the most stud¬ 
ied and utilized. Here, the interaction between the light 
and the free charges of the metal is exploited to achieve 
light confinement at a metallic surface and thus locally 
increase the electromagnetic field intensity [Ij. This ef¬ 
fect is used, for example, to increase the intensity of the 
Raman signal of a molecule or the efficiency of a fluo¬ 
rescent marker bonded to a target molecule 0,1- Yet, 
light confinement near the surface can also be exploited 
to control its propagation in the plane through a nano¬ 
structurization of the surface, which allows one to create 
waveguides and resonators with minimal fabrication re¬ 
quirements 0|. 

A promising strategy for such control and enhance¬ 
ment of light is based on photonic crystal (PhC) struc¬ 
tures supporting Bloch surface waves (BSWs). These 
are evanescent electromagnetic field modes that prop¬ 
agate along the interface between a periodic dielectric 
stack and an homogeneous medium [5|, Q. Their confine¬ 
ment relies on total internal reflection from the homoge¬ 
neous medium and reflection within a photonic gap from 
the multilayer. Although more complicated than metal¬ 
lic systems, PhC-based devices are highly customizable, 
and they do not suffer from the absorption losses that 
plague metallic structures M. 
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High optical quality multilayers can be grown by 
molecular beam epitaxy, but they are also commercially 
available from thin-film companies for many materials, 
from semiconductor to oxides. More recently, the gen¬ 
eral features of guided BSW have been characterized the¬ 
oretically and demonstrated experimentally IT2I ITi! . All 
these results seem to suggest the opportunity for a gen¬ 
eration of optical devices based on BSWs. This could 
be a very powerful approach for the development of in¬ 
tegrated optical sensors, but also a new and versatile 
platform to probe interaction between photonic modes 
and electronic nanostructures in classical and quantum 
optics, for instance in fundamental studies of the light- 
matter interaction in the strong-coupling regime |15l-|l9j . 
All of this requires the demonstration of optical res¬ 
onators, which are fundamental building blocks of any 
photonic platform, but are still lacking for BSWs struc¬ 
tures. 


In this work, we investigate the feasibility of realiz¬ 
ing a ring resonator on the surface of periodic dielec¬ 
tric stacks supporting BSWs. In this structure the ver¬ 
tical confinement is provided by exploiting the surface 
state, and the lateral confinement is defined by a nano¬ 
fabrication of the surface, which could be done using 
either standard lithographic approaches or low-cost im¬ 
printing techniques. In Ref. [13[ Liscidini showed that 
this structure is extremely flexible and can support a 
variety of guided modes depending on the structure pa¬ 
rameters and constituent materials. Very recently, Wu 
et al. reported on the fabrication and characterisation 
of bent waveguides supporting BSWs, demonstrating the 
feasibility of this geometry for more complex structures 


Here we shall consider TiO-zlS 1 O 2 multilayers and 
Polymcthyl-methacrylate (PMMA) ridges to work in the 
visible spectrum, for which the problem of designing 
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small and integrated resonators is typically more de¬ 
manding than in the infrared wavelength range due to 
the lack of transparent materials that can guarantee a 
large refractive index contrast. In this case one can work 
by etching a PhC crystal cavity in a membrane, but the 
resulting structure is usually very fragile [20, [2l[ ■ Alter¬ 
natively, there have been demonstrated Hydex or SiN 
high-quality ring resonators embedded in a silica matrix 
[22, [231 , but in this case the region in which one experi¬ 
ences the largest light-matter interaction is not directly 
accessible. 


While most of the theoretical tools necessary for the 
design of multilayers are already available, the design 
and optimization of BSW waveguides and ring res¬ 
onators is challenging. On the one hand, as the multilay¬ 
ered structure can be characterized by numerous layers, 
even a hundred, whose thicknesses can be as small as a 
few tens of nanometer, FDTD approaches can be largely 
time consuming and difficult to implement when the sim¬ 
ulation cell is particularly large and high spatial resolu¬ 
tion is required, as in this above-mentioned case. On the 
other hand, the typical size of the simulation in the plane 
of the multilayer, several micrometers, makes it difficult 
to describe the structure also in the reciprocal space by 
using Fourier Modal methods. Along the years, several 
strategies to solve for the confined modes in bent waveg¬ 
uides have been proposed: based on a numerical solution 
of Maxwell equations in terms of an eigenvalue problem 
24], expansion in Hankel wavefunctions in the cladding 
25j , conformal mapping of the bent waveguide [26[ , per¬ 
turbative approaches [27[, beam propagation methods 
[2^ | . finite element discretization [§9j], etc... 

Here, the different mechanisms at the base of the con¬ 
finement of light in the vertical direction and in the plane 
of the structure suggest to follow a strategy based mainly 
on effective index approaches, which reduce the dimen¬ 
sionality of the problem and thus are able to quickly 
explore several configurations in the parameter space. 
In particular, our goal is to understand what are the 
main parameters that limit the quality factor achievable 
in an ideal BSW-based ring resonator when scattering 
losses due to fabrication imperfection can be neglected. 
Beside demonstrating the feasibility of a ring resonator 
approach, these results will serve as a guide in the design 
and development of a complete BSW-based platform. 

The paper is structured as follows: in Sec. [2] we 
describe the structure under investigation and the ap¬ 
proach to the calculation of the mode quality factor of 
a PhC ridge ring resonator; in Sec. [3] we focus on the 
analysis of the intrinsic losses of BSW waveguide modes 
supported by a PhC ridge on a finite multilayer; in Sec. []] 
we deal with the bending losses associated with the typ¬ 
ical lateral confinement that can be obtained for modes 
in PhC ridges; in Sec. Owe outline a general strategy 
to design BSW-based ring resonators. Finally, in Sec. [G] 
we draw our conclusions. 



Fig. 1. (Color online) (a) Sketch of the PhC ridge ring res¬ 
onator. (b) Cross section of the corresponding PhC ridge. 


2. Structure and theoretical approach 

The structure we have in mind is a ring resonator of ra¬ 
dius R fabricated on the top of a multilayer supporting a 
BSW (see Fig. EH a)). The light confinement in the verti¬ 
cal direction is due to the surface mode, while the lateral 
confinement is given by the PhC ridge, whose section is 
shown in Fig. mb). This device could be fabricated 
using different kinds of dielectric materials, from semi¬ 
conductors to oxides, in principle to operate at a given 
working point in the whole spectral region, with a band¬ 
width strongly dependent on the geometrical parameters 
and material choice. The structure we have in mind is 
a good compromise between micro-disk resonators [30| . 
which are typically realized in high-refractive index con¬ 
trast platforms, and micro-pillar resonators [31 j. which 
are less demanding in terms of refractive index contrast, 
but they are also more difficult to fabricate. The features 
of our structure are expected to be intermediate with 
respect to the above-mentioned systems: quality factors 
larger than those typically observed in micro-pillars (but 
smaller than what observed in micro-disks) along with 
flexibility in terms of fabrication and materials. In the 
following we will focus on a specific example to provide 
a description of the main features of this resonator. 

The multilayer is periodic, with the unit cell composed 
of two layers with thicknesses c?Ti0 2 = 0.085 ^m and 
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Fig. 2. (Color online) Intensity profile in the yz -plane of 
a TE-polarized guided BSW calculated at Ao = 0.630 ym. 
The PhC is made of IV = 10 periods of a TiOz/SiOz unit 
cell (dTi0 2 = 0.085 fim and dsi0 2 — 0.128 pm) and the first 
TiOz layer is truncated with d^ Tst = 0.010 pm. The ridge 
material is PMMA with a rectangular section (0.800 pm x 
0.220 /im). The profile is obtained by means of a FDTD 
numerical simulation. 


dsio 2 = 0.128 /im, respectively. The multilayer has a fi¬ 
nite number N of periods and is truncated with the first 
layer made of TiOz with thickness dfi rs t = 0.010 /jm. 
The structure is designed to operate in the visible spec¬ 
tral range at Ao = 630 nm (about 1.97 eV) with riTi0 2 = 
2.58534 and nsio 2 = 1.54270. The choice of these ma¬ 
terials for the multilayer is convenient from a techno¬ 
logical point of view, for they are commonly employed 
in the fabrication of optical filters by sputtering, thus 
this plan structure is commercially available. Finally, 
we consider a PMMA (npMMA = 1.48914) ridge with 
height h = 0.220 pm and width w = 0.800 yon on the 
top of the multilayer. 

The structure shown in Fig. mb) is symmetric 
upon reflection with respect to the xz-p\ane. Conse¬ 
quently, the guided modes can be classified according 
to the eigenvalues of the mirror operator a xz in TE 
(Transverse-Electric, a xz = —1) and TM (Transverse- 
Magnetic, a xz = +1) modes. In the following we shall 
focus only on the fundamental TE mode, whose intensity 
profile, calculated by FDTD [H, [Hj], is plotted in Fig. 
[2j This mode is a perturbation of the BSW supported 
by the bare multilayer, and light is guided within the 
ridge likewise in a rib waveguide. It should be noticed 
that, unlike truly guided BSW, this mode can exist for 
any choice of the ridge height, with its pro perties and 
modal volume depending strongly on h [131 ] . Here the 
ridge height has been chosen to guarantee a good lat¬ 
eral confinement and, at the same time, keep the modal 
volume as small as possible. 

One can obtain a ring resonator by bending the 
straight waveguide onto itself in a circular shape (see 


Fig. Oja)). The resonant frequencies are obtained by 
the usual condition Q 

27T 

kL = —n mode 27rR = 2 tt m, (1) 

Ao 

where n mo de is the real part of the effective index of the 
mode in the bent waveguide, R is the ring radius, and 
m is an integer. In the absence of the x^-plane mir¬ 
ror symmetry, a rigorous TE/TM classification of the 
guided modes is no longer possible. However, in waveg¬ 
uides with a rectangular cross section, the inodes hav¬ 
ing the E r component dominant over E^ and E z (35| 
can be referred to as TE-like , and in the following we 
shall always consider these class of solutions. Finally, 
this structure is characterized by a continuous rotational 
symmetry and therefore, by choosing cylindrical coordi¬ 
nates (r,cj),z), the angular dependence of the fields can 
be written in the form e lm ^, with m integer. This prop¬ 
erty turns out to be extremely useful in the numerical 
description of the structure, as it allows one to reduce 
the dimensionality of the electromagnetic problem, for 
example in FDTD calculations. 

Our goal is the design of a ring resonator working 
at Ao and the computation of the mode quality factor, 
which can be used to measure the performances of the 
device. In the case of a BSW-based ring resonator, this 
task is all but simple, as an analytical approach is not 
feasible, and even a brute force 3D numerical calcula¬ 
tion, for example by using a FDTD method, can be 
extremely time-consuming on a standard personal com¬ 
puter or medium-sized server. Thus, in the view of opti¬ 
mizing the structure, which depends on a large number 
of parameters (unit cell composition, multilayer trunca¬ 
tion, height and width of the ridge, ring radius, etc...), 
one should rather look for an approximated but faster 
strategy. 

The structure and the mode profile in Fig. [2] seem 
to suggest a possible path towards the solution of the 
problem. Indeed, the strong field confinement of light in 
the ridge waveguide indicates that the resonant modes 
properties can be found by decomposing the 3D quest 
in two independent ID and 2D problems. First, we de¬ 
scribe the field confinement due to the surface state in an 
effective multilayer, then we find the resonant modes by 
modelling an effective 2D system using a FDTD method, 
which allows us to take into account also for the depen¬ 
dence of n mo de on the ring radius R. Similar approaches 
have been used to study lower-order !36j and whisper¬ 
ing gallery modes j3ll] in micro-pillar resonators, demon¬ 
strating surprising effectiveness even for R of a few pm. 
Very recently, the same strategy has been successfully 
applied to describe light confinement and propagation 
of BSWs in bent waveguides [lij . 

It should be noticed that the mode quality factor Q 
depends on the light leakage through the substrate, due 
to the finite number N of periods in the multilayer, and 
the bending losses, which, for a given ridge, are a func¬ 
tion of the ring radius R. One can start by considering 
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that Q measures the energy loss experienced by the elec¬ 
tromagnetic mode per cycle and, assuming that the two 
loss mechanisms are independent, can write 

1 _ 1 1 

Q ~ Ql + Q? (2) 

where Q± is the quality factor determined by the out-of- 
plane losses and Q m is the quality factor limited by the 
in-plane losses. Thus, first we deal with the description 
of light leakage in the substrate, then we compute Q\\. It 
should be noticed that here we are neglecting extrinsic 
losses, such as scattering at the interface between lay¬ 
ers due to surface roughness or fabrication imperfection. 
In a real device these contributions can be important, 
depending on the accuracy in the fabrication processes. 

3. Multilayer design and Q± 

Here we consider a sufficiently large ring radius R so 
that the mode field profile and propagation wavevector 
are slightly affected by the curvature, i.e. the effective 
index n mo de in © is close to that of a straight waveg¬ 
uide (see also Sec. G|), and the vertical and horizontal 
dynamic can be decoupled. For this reason, the quality 
factor Q± can be computed by determining the propaga¬ 
tion losses of the corresponding guided BSW (GBSW) 
in a straight ridge waveguide having the cross-section 
shown in Fig. [ljb). In particular, if one assumes a semi¬ 
infinite PhC, the mode is truly guided as there are no 
propagation losses. On the contrary, in a more realis¬ 
tic situation, when the multilayer has a finite number 
N of periods, the GBSW is characterized by a finite 
propagation length, even when scattering losses can be 
neglected, as light can tunnel trough the multilayer into 
the substrate. 

The analysis of light confinement in the yz plane for a 
straight waveguide could be performed by using FDTD 
calculations or other 2D numerical approaches. Yet, in 
the spirit of suggesting a reliable but also fast approach 
that would allow one to rapidly explore the parame¬ 
ter space of the structure, we adopt an effective index 
method (EIM) by separating the field confinement in 
the y and z direction. It should be noticed that this 
approximation has already been proved to be reliable to 
study analogous structures [HI, l37M40| . More recently, 
the accuracy of this strategy for the description of guided 
modes in PhC ridges has been investigated in Ref. (l3 |. 
but the issue of propagation losses has not been dis¬ 
cussed. 

The EIM is conveniently employed to decouple the 
field confinement in 2D electromagnetic problems, pro¬ 
vided that light is well confined in the ridge at least in 
one dimension. Thus a 2D problem is reduced to the 
composition of two ID tasks, easily solved by means of 
the transfer matrix method (TMM). 

First, the EIM fictitiously collapses the guiding 2D 
structure in the direction where the most intense light 
confinement takes place. This is obtained by dividing 
the structure in ID multilayers (in accordance with the 


guiding regions shape) and calculating in each region the 
confined mode profiles and their effective indices. Subse¬ 
quently, the structure in each slice is replaced by an ho¬ 
mogeneous layer with a refractive index n e ff correspond¬ 
ing to that of a proper confined mode. Therefore, one 
can treat the most significant light confinement mecha¬ 
nism in each section exactly, and obtain an effective ID 
structure along the complementary direction, where the 
less intense localization of the field takes place. This al¬ 
lows the computation of the 2D field profile, the mode 
dispersion relation, and the propagation losses. 

One should note that the EIM predictions, although 
approximated, are generally very close to the values cal¬ 
culated by more time-consuming 2D numerical methods. 
However, the goodness of this approach depends on the 
structure geometry and field distribution. 

Referring to Fig. [2] it is clear that in our case the 
field confinement is strong in both y and z direction. 
For this reason, we expect that dividing the structure in 
the vertical or in the horizontal direction (although this 
implies different approximations) should produce quite 
similar and accurate results. In the following we will 
study the properties of the confined modes supported 
by our structure employing both the strategies. 

3.A. Horizontal EIM 

Here we exploit the strength of the vertical confinement 
of the field, due to the PBG mechanism on the PhC 
side and the TIR phenomenon at the P MM A/ Air in¬ 
terface, and hence we divide the structure along the z 
direction. Following the EIM described above, one sub¬ 
stitutes the PMMA ridge with an effective homogeneous 
layer of thickness h and refractive index n e ff, s iab(w), cor¬ 
responding to that of the fundamental TM guided mode 
supported by a symmetric PMMA slab waveguide of 
width w in air. Note that here TM refers to the plane 
of the PMMA slab waveguide. This is equivalent to the 
horizontal cross section of the PMMA ridge (see Fig. 
[3))a)), thus we shall refer to this method as the hori¬ 
zontal EIM (HEIM). The effective index can be writ¬ 
ten as n e ff,siab = Keff,siabc/w, where K e ff,siab is the zero 
of the T 22 ,siab(fc, w) element of the slab transfer ma¬ 
trix. Finally, one can calculate the dispersion relation 
of the GBSW by finding the guided modes of this effec¬ 
tive multilayer, searching for the solutions of the equa¬ 
tion T 22 ,mui(fc, w) = 0. In general, these are solutions 
of the kind kgsbw = /3 + ry, where (3 = (w/c)riGBSW 
and 7 = 1/ (2L pmp ), with T pr0 p the mode propagation 
length. 

In Fig. fflb) we show the real part of the dispersion 
curve corresponding to our parameters. As expected, 
the mode dispersion relation is within the PBG and be¬ 
low the air light line. In particular, for Ao = 630 nm 
we have (3 = 13.5164 /zm _1 and ugbsw = 1-3553; the 
corresponding mode intensity profile is plotted in Fig. 
He). 

For a sufficiently large N (typically more than 5, de¬ 
pending on the refractive index contrast of the materi- 
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Fig. 3. (Color online) (a) Sketch of the dimensionality re¬ 
duction of the electromagnetic problem by using the HEIM, 
(b) Real part of the GBSW dispersion relation. The yellow- 
shaded regions represent the PBGs. The dispersion relation 
is calculated considering the modal dispersion in the first 
effective layer, (c) E y held intensity in the truncated multi¬ 
layer. 


als), the dispersion relation is independent of the num¬ 
ber of periods of the multilayer. On the contrary, the 
mode propagation length depends exponentially on TV, 
for Lp r0 p is directly related to the imaginary part a of 
the Bloch wavevector, which controls the exponentially 
damping of the electric field in the multilayer. In gen¬ 
eral, one has 

iprop oc i exp[2TVaA], (3) 
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Fig. 4. (Color online) Propagation length L as a function 
of the number of periods N in the multilayer, computed by 
means of the HEIM. On the right axis the corresponding per¬ 
pendicular quality factor Q± (when Fig. [U(b) is the section 
of a ring resonator) is reported. 
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also write Q± = L prop /X = P/(Att^/). In Fig. 0]we show 
Q± calculated from the GBSW complex wave vector. As 
expected, we have an exponential growing of Q± with TV. 
It should be noticed that in our structure one can ob¬ 
tain propagation distances of about 10 m with only 10 
periods. In particular, this increasing is well fitted by 
a curve of slope (2aA) as suggested by Eq. (J3]). Thus, 
our results show that a fast and reliable optimization of 
the multilayer can be done very easily by computing the 
imaginary part of the Bloch wave vector as a function of 
the sole unit cell composition. 

In Fig. [5] we plot the imaginary part of the Bloch 
wavevector as a function of /3 at Ao = 630 nrn for 
our structure parameters. We also indicate the work¬ 
ing point corresponding to our final choice of the multi¬ 
layer truncation and the ridge parameters. It should be 
noticed that this point does not coincide with the curve 
maximum, which would guarantee the fastest held decay 
in the multilayer. Indeed, although one usually operates 
in the region of large a , the BSW ring resonator is de¬ 
signed to minimize also the bending losses in the ring, 
which of course depend on the ridge parameters. Thus, 
in general, the optimal parameters do not necessarily 
lead to the fastest growing of Q± with TV, due to a com¬ 
promise between vertical and lateral light confinement. 


where a is 


a = — arccosh 
A 




(4) 


in which one has applied the Bloch-Floquet theorem to a 
generic PhC with period A and transfer matrix T\ [3T| . 

Since the quality factor of the resonator can also be 
viewed as the number of held oscillations in the ring be¬ 
fore the held intensity is reduced by a factor 1/e, we can 


3.B. Vertical EIM 

The held intensity prohle reported in Fig. [2] shows 
that a strong conhnement of the GBSW takes place 
also in the y direction, due to the TIR mechanism at 
the PMMA/Air interfaces. Thus, one can obtain re¬ 
sults similar to those shown in the previous section fol¬ 
lowing the EIM prescription and dividing the structure 
vertically in two different regions, named bare , where 
the multilayer cladding is simply air, and loaded , in 
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Fig. 5. (Color online) Imaginary part of the Bloch wavevec- 
tor in a PhC with T 1 O 2 ISiOi unit cell (dTio 2 = 0.085 /im, 
d Si o 2 = 0.128 fim), as a function of the wavevector. The 
curve is calculated at Ao = 0.630 /xm. The working point 
corresponding to our choice of the multilayer truncation and 
the ridge geometry, calculated both with the HEIM and the 
VEIM, is reported. 


which the multilayer has an additional layer of height 
h corresponding to ridge region (see Fig. [Ha)). We 
search for the TE surface states supported by each mul¬ 
tilayer, and we calculate the corresponding effective in¬ 
dices n e ff,bare and loaded at A 0 . The straight waveg¬ 
uide is finally approximated by a vertical slab of index 
n-eH ,loaded with cladding of index n e ff jb are- Finally, one 
searches for the fundamental TM mode, whose wavevec¬ 
tor is kgbsw = /? + * 7 - Note that here TM refers to the 
plane of the effective slab. We shall refer to this second 
effective index approach as the vertical EIM (VEIM, see 
also Ref. [HI). 

The dispersion curves associated to the bare and 
loaded multilayers are reported in Fig. ( 6 l)b) , and 
the effective indices calculated at Ao = 630 nm for 
N = 10 are n e ff,bare = 1.1225 + 2.0602 ■ 10 -7 * and 
n e s, loaded = 1.3853 + 2.2347 • 10 -9 *, respectively. The re¬ 
sulting GBSW wavevector is kgbsw = (13.4332+1.2205- 
10 ~ 7 ) /xm _1 . 

In Fig. 0 we show Q± calculated by means of the 
VEIM as a function of the number of periods in the PhC. 
Similarly to the HEIM case, Q± increases exponentially 
with N , and the slope is very close to that shown in Fig. 
[H For a direct comparison, we report the working point 
obtained with the VEIM in Fig. E] along with the HEIM 
one. The two methods show a good agreement in the 
prediction of both the propagation constant [3 and the 
losses of the GBSW. 

4. Ring design and Qy 

Here, the calculation of the in-plane quality factor Qy 
is performed only by means of a combination of VEIM 
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Fig. 6. (Color online) (a) Sketch of the dimensionality reduc¬ 
tion of the electromagnetic problem by using the HEIM, (b) 
Real part of the dispersion curves for the bare and loaded 
modes. The yellow-shaded region represents the PBG. 
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Fig. 7. (Color online) Propagation length I as a function 
of the number of periods N in the multilayer, computed by 
means of the VEIM. On the right axis the corresponding per¬ 
pendicular quality factor Q j_ (when Fig. [TJb) is the section 
of a ring resonator) is reported. 
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and FDTD simulation in cylindrical coordinates, which 
is analogous to the approach adopted in Ref. [3| , where 
it is experimentally validated. First, the structure is di¬ 
vided in bare and loaded multilayer regions (see Fig. 
EKb)), with refractive indices n e ff,bare and n e ff,loaded, re¬ 
spectively. This allows us to approximate the 3D struc¬ 
ture with a 2D ring resonator (see Fig. [51(a)). 





Fig. 8. (Color online) Top view of the circularly bent effec¬ 
tive slab waveguide. The light blue region refractive index is 
u e ff, loaded = 1.3853 while the surrounding medium is charac¬ 
terized by U e ff bare = 1.1225. 

The computation of Q m around Ao is performed using 
Meep, a free FDTD simulation software package devel¬ 
oped at MIT j4l|. This code is characterised by the 
possibility of exploiting the continuous rotational sym¬ 
metry of the ring by working in cylindrical coordinates 
and solving the angular dependence of the fields analyt¬ 
ically, according to e* m< ^, with m integer. This reduces 
the calculation to a simple ID simulation, with large 
benefits in terms of the computational effort. 

This approach takes into account the modification of 
the mode effective index due to the waveguide bending. 
It is interesting to compare the results calculated for 
the straight ridge by means of both the HEIM or VEIM 
with the value of the mode effective index in the case 
of a bent waveguide. The comparison is shown in Fig. 
El As expected, for larger R, the effective FDTD results 
approach the mode index calculated with the VEIM. It 
should be noticed that in the range of R considered in 
our work all the three values are very close, with differ¬ 
ences corresponding to an error in the determination of 
the resonance frequency less than the free spectral range 
of the ring resonator. This confirms the goodness of our 
approximation for the parameters under consideration 
and highlights the possibility of determining the reso¬ 
nance frequency (but not the losses) by means of the 
analysis of the sole straight waveguide, even in the case 
of relatively small R. 

In Fig. [10] we show Qy as a function of the ring ra¬ 
dius. As expected from a standard treatment of bending 
losses, Q\\ increases about exponentially with R. We also 
show Q± versus N, so that one can immediately see the 
essential requirements in order to meet a given Q target. 

It should be noticed that the value of the quality fac¬ 



Fig. 9. (Color online) Effective index of the GBSW calculated 
by means of the VEIM and 2D FDTD as a function of the 
ring radius R. The values calculated with the HEIM and 
the VEIM for the corresponding straight waveguide are also 
shown. 
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Fig. 10. (Color online) Lateral quality factor Qy as a function 
of the bending radius R calculated by means of the VEIM and 
FDTD numerical simulations. The result is compared with 
the perpendicular quality factors Q± obtained in section [3] 

tors in such a small ring resonator exceeds those that are 
experimentally observed in structures characterized by 
a stronger refractive index contrast [42| or larger radius 
R 0. Naturally, these theoretical values are obtained 
by neglecting scattering losses and fabrication imper¬ 
fections, which reduce the quality factor dramatically. 
However, the quality factors calculated here are so high, 
that even if scattering losses (and every other extrinsic 
loss mechanism) will affect the resonator considerably, 
the actual device would still display good light confine¬ 
ment properties. For instance, let us consider a ring 






















resonator with N = 10 periods and radius R = 9 /. im. 
For this structure the quality factor is Q ~ 10 7 and 
even allowing the scattering losses to lower this value by 
three orders of magnitude, we would still obtain a ring 
resonator with Q ~ 10 4 within a footprint of ~ 250 pm 2 . 

5. Design strategy and final remarks 

The design of a BSW-based ring resonator to achieve a 
target quality factor is critically dependent on the con¬ 
stituent materials, as well as on the PhC truncation, 
and on the ridge section. For this reason, it is useful 
to outline a basic strategy for designing and optimizing 
BSW-based resonators. 

In the previous sections we have considered a specific 
Ti0 2 /Si0 2 BSW ring resonator working at A 0 = 630 
nm. Although the final structure depends critically on 
the constituent materials and target wavelength, it is 
possible to outline a general design strategy to obtain 
structures to work in a different spectral region and char¬ 
acterised by a target quality factor and modal volume. 
The latter is essentially determined by the exponential 
decay of the electromagnetic field in the multilayer and 
the ring radius. Despite the main goal of this work is 
not the optimization of a structure for a specific applica¬ 
tion, one can be guided in the design of new BSW ring 
resonators by a few simple rules. 

As we have seen, a BSW ring resonator is made up of 
several layers and different materials, and it is character¬ 
ized by a rather complex geometry that makes the struc¬ 
ture modelling quite challenging. Indeed, it requires to 
perform several numerical simulations (looped over all 
the variables) to cover the vast parameters space. In 
this regard, the EIM-based approach we have proposed 
here is actually a valuable tool, suitable for this pur¬ 
pose. It allows to obtain the relevant physical properties 
of the confined electromagnetic modes introducing min¬ 
imal errors while requiring a very small computational 
cost. Yet, it should be noticed that, even in this case, 
one cannot explore such a large parameter space at the 
same time. 

Once the target wavelength is chosen, one starts by 
optimising the unit cell of the multilayer that will sup¬ 
port the entire structure. Given the materials compat¬ 
ible with the most accessible fabrication technique and 
the specific application one has in mind, the best strat¬ 
egy is looking for strong refractive index contrast and 
transparency. Indeed, this will lead to a large PBG that 
offers flexibility in the design of the PhC ridge and res¬ 
onator. This will also determine high Qj_s with small N, 
for the larger is the index mismatch in the unit cell, the 
faster is the field exponential decay in the multilayer. 
For the choice of the layers’ thickness, a good thumb- 
rule is considering quarter-wave stack, where this condi¬ 
tion has to be satisfied at the working point below the 
light line corresponding to the upper cladding material 
and not at normal incidence as it is usually done in dis¬ 
tributed Bragg reflectors. 

The second step is choosing the multilayer truncation 
and ridge parameters , which determine the main prop¬ 


erties of the mode supported by the PhC ridge. In this 
communication we have dealt with a mode characterised 
by a strong field confinement in the ridge region. How¬ 
ever, a PhC ridge can support a larger variety of modes, 
which all could be used for the realization of BSW ring 
resonators (see e.g. 0 )- 

With respect to the case investigated in this paper, 
the position of the bare and loaded modes inside the 
PBG at a given Ao is crucial in the optimisation of 
both Q± and Qy. On the one hand, the best vertical 
confinement of the surface state is obtained when both 
modes are close to the PBG center (where a is highest). 
On the other hand, one looks also for a large contrast 
A n eS = n e ff,loaded - n e ff,bare, which guarantees better 
lateral confinement and small bending radii. These two 
conditions point in opposite directions in the parameter 
space, thus one looks for a compromise between lateral 
and vertical confinement. In our case, without any par¬ 
ticular application in mind, we have designed our struc¬ 
ture to maximize the effective index mismatch by mak¬ 
ing the bare mode effective index close to the air light 
line; the loaded effective index is engineered to be as 
close as possible to the PMMA light line while main¬ 
taining at the same time the vertical waveguide single 
mode. This results in a very compact BSW ring res¬ 
onator, characterised by large quality factors. 

The last step is the structure trimming by using 
a more accurate fully-3D simulation, for example by 
means of FDTD. Yet, depending on the size of the fi¬ 
nal structure, this might be hard, if not impossible, on 
a standard personal computer. 

6. Conclusions 

In this work we proposed an original resonator consist¬ 
ing in a photonic crystal ridge supporting a guided BSW 
and bent in a circular shape. We have studied its opti¬ 
cal properties, including the dependence of the quality 
factor on both the ring size and multilayer composition. 
We outlined a very general strategy that could be used 
to design and optimize BSW ring resonators made of 
different materials and operating in different spectral 
regions. This approach is based on assuming vertical 
and in-plane optical losses as independent, and comput¬ 
ing the corresponding perpendicular quality factor Q± 
and in-plane quality factor Qy by means of TMM and 
effective 2D FDTD calculations. As an example, we de¬ 
signed and analyzed a BSW ring resonator constituted 
of a PMMA ridge waveguide on a truncated periodic 
Ti0 2 /Si0 2 multilayer, operating at Ao ~ 0.630 /im, 
having a theoretical value of Q exceeding 10 7 , and a 
device footprint of ~ 250 pm 2 . 

Given the strong interest in BSWs for the realization 
of integrated optical sensors and as a tool to investigate 
the light-matter interaction at the fundamental level, 
and the recent experimental demonstration of guided 
BSWs in bent waveguides, we believe these results might 
be of great interest to that part of the community work¬ 
ing in integrated optics. 
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